Vector vibrations and the loffe-Regel crossover in disordered lattices 



S. N. Taraskin and S. R. Elliott 
Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CBS lEW, UK 

(Dated: February 1, 2008) 

The spectral density for vector vibrations in the f.c.c. lattice with force-constant disorder is 
analysed within the coherent potential approximation. The phase diagram showing the weak- and 
strong-scattering regimes is presented and compared with that for electrons. The weak-scattering 
regime for external long-wavelength vibrational plane waves is shown to be due to sum-rule corre- 
lations in the dynamical matrix. A secondary peak below the Brillouin peak for sufficiently large 
wavevectors is found for the lattice models. The results obtained are supported by precise numerical 
solutions. 



I. INTRODUCTION 



Vibrational excitations in disordered structures is an active field of investigation by different metJaftds. The first 
approaches were developed many decades ago for the case of substitutional alloys (see e.g. Ref.UBBcl and refer- 
ences therein). Disorder was introduced there by the random positioning of two or more atomic snecies (of different 
masses) onto lattice sites. Successful analytical methods have also been adapted for such problcmsDElfl. Tjiie rigidity- 




percolation problem has been solved to some extent numerically and within-tiLe mean-field approximatiopi^prjl. y^tpmic 
vibrations in [Stpti&tures with fractal geometries were investigated in Rcfs.ElEJ. Some phenomenologi( 
more realistidlSOO models of atomic vibrations in glasses have been developed in order to describe the thermodynam- 
ical and phonon-transport peculi|arities of such materials. Various computer experiments have also been performed 
on disordered structures (see e.g.ES and references therein). 

Recently, the methods of the mode-coupliag theory have been applied to investigate short-time vibrational dy- 
namics, especially in the low-energy regimellj. The low-energy behaviour of the vibrational density of states and 
dynamical structuije factor have been investigated in disordered lattice models within the coherent potential approxi- 
matipi|ir-(f^|4|)ElL2lc3, and in positionally disordered models treated perturbatively in the low inverse-particle-density 
limitE3'MPlo. The spectral properties of the Euclidian random matrices associated |jsdth topologically disordered 
systems^ have been investigated within a^upersymmetric statistic fieid theory in Ref.EIl. The spectral properties of 
the Laplacian defined on a random graphca and " small- world" latticescj are related to topologically disordered solids 
and have also been investigated. 

Previously, vibrational displacements in disordered structures have been treated as scalars, which is a convenient 
simplification of the|j/ector nature of atomic vibrations (although vector vibrations were considered in the rigidity- 
percolation problemQ). However, a consistent treatment of vector vibrations in disordered lattices has not been 
performed up to now. The main purpose of this paper is to fill this gap and to give a mean-field description of vector 
vibrational models, namely: (i) to derive the CPA equations for vector vibrations in disordered lattices with force- 
constant disorder; (ii) to solve them analytically in limiting cases and numerically through the whole energy range; 
(iii) to apply the solution to investigate the scattering properties of external plane-wave vibrational excitations and, 
in particular, the loffe-Regel crossover from a weak- to strong-scattering regime; (iv) to demonstrate the importance 
of sum-rule correlations in the dynamical matrix for scattering properties of vibrational external plane waves. One of 
the main results of the paper is the phase diagram for vector vibrations in the representative f.c.c. lattice showing the 
ranges of weak- and strong-scattering regimes for different degrees of disorder. Another result concerns a comparative 
analysis of the vibrational problem with the electronic one, and an explanation of the essential differences in phase 
diagrams for these problems. All the major derivations and conclusions are supported by precise numerical solutions 
for the same model, which demonstrate the very good performance of the mean-field CPA approach, both for the 
density of states and the spectral density in a plane-wave basis (dynamical structure factor) in the whole energy range, 
except for the high-frequency tail region containing localized states. 

The rest of the paper is arranged in the following manner. The model is defined in Sec. IL The properties of the 
effective mean field are analysed in Sec. IIP The spectral densities and the scattering properties of the external plane 
waves are discussed in Sec. IV. The conclusions are given in Sec. 0. 
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II. MODEL 



A. Hamiltonian 



Atomic vibrations in the harmonic approximation can be treated in the classical limitu. Thus the problem is 
reduced to the Hamiltonian formalism, in which the dynamical operator plays the role of the Hamiltonian, and the 
squared frequency, w^, is an energy, e = uP' . The relevant Hamiltonian (dynamical operator) for vector vibrations of 
'N interacting atoms of unit masses (m^ = 1) can be written in the following form: 



H 



(1) 



where we have used the site orthonormal basis of vectors |?,a), the only non-zero component of which is related to 
the displacement of atom i (i = 1, . . . , A'^) along the Cartesian direction a (a = 1, ... , D). The (dynamical) niatrix 
elements Hiaja' , being the second partial derivatives of the interatomic potential, obey the following sum rule£]: 



Hi, 



(2) 



These sum-rule correlations between the elements of the D x D diagonal and off-xliagonal blo cks d istinguish the 
vibrational problems from the electron problems described by similar HamiltoniansoQ (see Sec. IV E). The matrix 

elements of the Hamiltonian (|^), generally speaking, are functions of all equilibrium atomic positional Vjectors, rf^\ 

Our aim is to solve the eigenproblem for the Hamiltonian (|l|). This can be easily done for a crystal^ but is not 
so simple for a disordered structure. The analysis starts with a choice of the vibrational model for a disordered 
system. This could be: (i) a crystalline, latt i ce , o f atoms of random masses (mass disorder) and/or connected by 
random springs (force-constant disorder )EjG2I'EHIe1I; (iil-xinaitionally disordered structures created, for example, by 
random atomic displacements around crystalline sitesoOEj, or by other ways, e.g. the Bethe lattice or by bond 
switching^; (iii) atomic structures, created e.g. by molecular Hyngrnirfj; frlat;j;ira1 or ab-initio), in which intrinsic 
positional disorder defines the force-constant disorder (see e.g. Refs.E3o'L£(L3E3). The simplest models are from the 
first class. Analytical approaches can be easily developecL£ox_|them due to the existence of a well-defined reference 
crystalline system (i.e. the same lattice without disorder )o'B'LrL3. Surprisingly at first sight (see below), lattice models 
with force-constant disorder mimic very well the main features of topologically disordered glasses and even liquids, 
e.g. the boson peak, l ^|lj| . e^ j j^*^ tails, main peaks in the density of states and zero-energy singularity in the 



instantaneous spectrumEiEilt3cJ'c2l. The know^Ajawback of lattice models is related to the fc^-frequency dependence 
of the width of the dynamical structure factorESEilEa'cll as compared to the parabolic behaviour found experimentally 
(see e.g.cj). This fc^-dependence leads to a parabolic frequency scaling of the dynamic structure factor, S{k,uj\ 
for cj — > 0, which is in contrast to the frequency-independent S{k,uj) found for topologically disordered modelst^lt 



B. Distribution of force constants 



In this paper, we consider as a vibrational model the crystalline f.c.c. lattice of atoms, in which the nearest 
neighbours are connected by springs that are unstretched in equilibrium. The matrix elements of the Hamiltonian (|l|) 
in this case have the following form: 

Hiaja' — ^'^iji^ij'^)ai^ij'^)a' : if * 7^ J i (3) 

where f,-^^ — (r'"'' — r ■'^^) /Ir^"-' — r,-"'] is the unit connection vector between atoms i and j in equilibrimn. The 
spring constants, Kij, are supposed to be independent of atomic positions and are random values taken from a 
certain probability distributio n , i p(jf) ti D ifferent functional forms of p{k) have been used for various problems, e.gp-a 
cornbination of two (5-p£unctionf£IleJc3'c£l, (also used in the rigidity-percolation problemEl) and the normal (gaussian)E3, 
boxLJ and hyperbolicE^ distributions for glassy problems. Below we use the box distribution: 

p(«) - ^ [0{^ -Ko + A)- e{K -Ko- A)] , (4) 

where the Heavyside function 9{x < 0) = and 9{x > 0) = 1, kq is the spring constant in the reference crystal and 
2A is the full width of the distribution. The particular form of the peak-shaped distribution is not important for all 
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results presented below from a qualitative point of view, as we have checked by using the normal distribution instead 
of the box one. The box distribution also has the advantage of not needing a lower cut-off to ensure positive values 
of force constants, as in the case for the normal and hyperbolic distributions. 

If the distribution given by Eq. (^) is wide enough, i.e. A > kq (kq = 1 below), then negative spring constants 
appear in the system which could result in vibrational instability of the model in equilibrium, i.e. the occucxepce of 
negative eigenvalues. Such a situation is of direct relevancec3 to an instantaneous mode analysis in liquidalIrc3CJ. 

The vibrational spectrum of a disordered system is defined by the properties of the Hamiltonian (dynamical) 
matrix. In the case of a crystal, the Hamiltonian matrix is well structured (e.g. cyclic band diagonal). For topologically 
disordered glasses with conserved local order, the Hamiltonian matrix is structured similarly to that for the crystalline 
counterpart (with the same local order) but now, due to the presence of topological disorder: (i) the positions of 
elements are distributed about their positions in the crystalline matrix (locational disorder) and (ii) the values of the 
elements are also distributed because of different interatomic distances. The first effect is completely ignored in the 
disordered-lattice models and the total influence of disorder is attributed entirely to the second effect. The question 
is, how good is this approximation? 

One criterion relates to similarities in the probability distributions for the values of the matrix elements for both 
types of models. We have compared the distributions of the Hamiltonian matrix elements for two vibrational models: 
(i) the f.c.c. disordered-lattice model discussed above and described by the matrix elements HI^^^, obeying Eq. (|^) 
and (ii) a topologically disordered model of a single^oinponent glass with predominantly icosahedral order (IC-glass) 
obtained by classical molecular dynamics (see Refs.EJ'EB for more detail). In the latter model, the Hamiltonian matrix 

elements, Hf^^l^,, are defined by the pairwise interatomic potentially, T^(|r|^^|), its spatial derivatives, V^'(|r|?^|) and 



V" 



and the relative atomic positions. 



rr glass 

ia,ja' 



'i^ii )a' 



'1 ^' 1^(0), 



V'i\r[f\) 
1.(0) I 



(5) 



The Hamiltonian matrix contains D x D diagonal and off-diagonal blocks, the diagonal {Hia^ia and Hia,ja) and 
off-diagonal {Hia,ia' and Hia,ja') elements being distinctively distributed. For the diagonal elements in the off- 
diagonal blocks, it is convenient to compare the traces taken with opposite signs. They represent the force constants, 
4* = Ea^lajc and 4'"^^ = Ec^fa?," = {y"{\^f\) + 21/'(|rg^|)/|rg)|), for interactions beween atoms i and j in 
lattice and glassy models, respectively. Similarly the traces of diagonal blocks correspond to the total force constants, 
K^f* and nfi'^^'', for atom i in lattice and glassy models, respectively, where ^'^^feiass) _ ^iat(giass) 

The four distributions for both models presented in Fig. |l|(a)-(d) show qualitatively the same features. The distri- 
butions of the diagonal elements in the diagonal blocks, pinfl^^^) and /9(K-f*), are both peak-shaped with remarkably 
comparable peak positions and peak widths (see Fig. |^(a)). The distributions of the off-diagonal elements in the diago- 
nal blocks are both centred around zero and have comparable widths (see Fig. |l|(b)). The distributions of the elements 
for the off-diagonal blocks are much more model dependent. The solid line in Fig. ^c) represents the distribution 
Pi'^fj^^'') all interacting atomic pairs. The distribution function has a complicated form, the features of which can 
be understood in terms of the functional forms of the interatomic potential and the atomic pair correlation function. 
For example, if nearest-neighbour interactions only are taken into account (as in the disordered-lattice model), then 
the distribution function p{Hfj^''^) has a single peak (see the inset in Fig. |l|(c)), which corresponds to the first peak 
in the pair correlation function^. The simplest model distribution which can mimic the function p{K^j^^) (at least 
its scale) is the box distribution (the dashed line in Fig. |l|(c)). The distribution of the off-diagonal elements in the 
off-diagonal blocks (see Fig. |l|(d)) has a maximum at zero and is symmetric about it. The disordered-lattice models 
are intrinsically anisotropic, so that the distribution function, p{H^^-^,), can look quite different to piHf^^'^^,). For 
example, for the ordered f.c.c. lattice (A = 0), it consists of two (5-functions symmetrically located around zero. With 
increasing disorder, these (^-functions are transformed into two box distributions which merge at A = 1, forming the 
single box distribution shown in Fig. |l|(d) by the dashed line. For larger A > 1, these two boxes overlap around zero, 
thus resulting in a box-shaped peak in this region, i.e. the distribution function piHj^j^,) becomes closer in shape 

to the peak-shaped function p{Hf^l^,). 

The above comparison demonstrates that the very simple disordered-lattice model with a box distribution of force 
constants is able to reproduce correctly the ranges of the distributions of the Hamiltonian matrix elements for a 
topologically disordered glass with minimal short-range order (i.e. the IC-glass). This leads to the plausible conjecture 
that the physical origins of different phenomena, such as the boson peak, loffe-Regel crossover and secondary peak in 
the dynamic structure factor (analysed below), localization-delocalization transition, etc., are common for topologically 
disordered glasses and disordered lattices. 
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The other important point to emerge concerns the possibihty of estimating the degree of disorder in disordered- 
lattice models at which they are able to mimic glassy models. For small disorder, all the distributions for the 
disordered-lattice model presented in Fig. || are quite narrow and are not comparable in width with the distributions 
for the glassy model. They become similar only for sufficiently large disorder (A ^ 1). 

Below we use the simplest (box) distribution of force constants which allows us to analyse the problem analytically 
and which mimics some of the features of topologically disordered glasses by disordered lattices. A natural extension of 
the lattice models is in the use of more realistic force-constant distributions for crystalline counterparts of topologically 
disordered glasses and a comparison of the resulting properties of such models with the corresponding glasses. A 
possible way to mimic disorder in glasses is to introduce thermal (positional) disorder in the_crystalline counterpart. 
We have done this for the IC-glass and its crystalline counterpart, the Frank-Kasper cr-phaseo. A cr-phase lattice has 
been heated by means of molecular dynamics and the dynamical matrix for the instantaneous cofiguration has been 
calculated. The elements of such a matrix do not exhibit locational disorder and only their values are distributed. 
Using temperature as a control parameter it is possible to obtain the distributions of matrix elements very similar to 
those ones for the topologically disordered IC-glass (see Fig. 6 in Ref.ES). The vibrational density of states (VDOS) 
for instfMijljMcous configurations of the cr-phase, as a consequence of such a similarity, resembles the VDOS for the 
IC-glassc3E3. This supports the use of lattice models for mimicking topologically disordered structures. We should 
mention, however, that the elements of the dynamical matrix (belonging to the Euclidian ensemblecj) for positionally 
disordered crystalline counterparts are correlated with each other within a certain row. These correlations are absent 
for the simple box distributions used below. 



C. Coherent potential approximation 



The solution of the eigenproblem for the Hamiltonian (|i|) is known for a crystal, when A = in Eq. (^). The 
eigenmodes then are plane waves |k, (3), characterized by the wavevector k, the branch number /? and the eigenenergies 
ek/3- The VDOS, g"^'^^{e) = {DN)~^ Sk/3'^(s ~ ^k/?), forms a band (see the dashed line in Fig. ||) containing van 
Hove singularities, and behaves as (e)'^^/^^~^ around the band edges0. Disorder (A > 0) brings new features to the 
spectrum (see Fig. (i) the van Hove singularities are washed out; (ii) a high-frequency tail containing localized 
states occurs; (iii) extra states appear in the low-energy regime. 

The exact solution of the vibrational problem with random spring const a^oi-s |is .not known, but a number of ap- 
proximate methods, both for electrons and vibrations, have been developedQ'ElE3C3. We are interested in the global 
features of the spectrum, except the tail region where the states are localizetLjAn adequate approach in this case is 
a mean-field treatment within the coherent potential approximation (CPA)o'Elt3, which we use below. The CPA is 
known to be very successful in a description of the electronic and vibrational properties of substitutional alloys and is 
expected to work very well in the whole energy range, except the far tails containing strongly localized states, where 
field-theoretical approaches are much betterEil. 

The Hamiltonian (^ can be rewritten in the bond representation. 



H 



(6) 



with \ij 

of disordered systems, g{e), 



a) — \j,a)), which is a convenient form for a subsequent mean-field treatment. The VDOS 



g(e) = — TrA(e) = 
' ND ^ ' 




\{d\k,(3)f5{e~ed) 



(7) 



((5(e — H)), where the angular brackets, (...), dene 
i|k, /3)p, in Eq. (|^) have a simple physical meaningt 



can be defined via the spectral-density operatorE^I, A(£) = 
averaging over random spring constants. The coefficients, 

they define the weight of the disordered eigenstate \d) in the crystalline state |k, f3), which is a plane wave of a certain 
polarization. The importance of the spectral-density operator in the plane-wave basis iSpUelated to the fa ct that the 
matrix elements Ai(;/3(e) = (k,5|A|k/3) are proportional to the dynamical structure factortj (see Sec. IVD) which can 
be measured experimentallyEJ. 

Both the VDOS, g{e), and the spectral density, Akfj{e), can be found for vector vibrations within the mean-field 
approximation (see Appendix A): 



g{e) = --Im 

TT 



1 



^cryst 



S 



(8) 
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and 



^k/3(e) 



'7r[e-z'(£)ek/3]' + [z"(£)ek/3]' 



(9) 



via the crystalline Green's function, G"^^^, and the complex effective mean field, z{e) = z'{e) + iz"{e) — K{e)/Ko (the 
dimensionless effective force constant), which can be found from the solution of the self-consistent Eq. (BqI). 



III. EFFECTIVE FIELD 

As seen from Eqs. (|8|)-(||), the spectral properties of a disordered lattice, within the CPA, are characterized by 
the energy dependence of the effective dimensionless force constant (effective field), z{e). In this section, we discuss 
the properties of the effective field obtained by the solution of the self-consistent equation (|2^) (see Appendix A) 
for different values of disorder and in different energy ranges. The most important regime from the viewpoint of 
comparison with experimental data for the dynamical structure factor is the low-energy limit (e — > 0), analyzed 
below. 

For the particular choice of the box-like force-constant probability distribution given by Eq. (^), the integration in 
Eq. (ES) can be done analytically, resulting in the following expression: 



1 



1 



2Aa(e 



dn 



1 



z(e) 



z(e) 



, 



(10) 



with 



a{e) = Ko{ij\G\ij) 



(11) 



Eq. ( |10[ ) should be solved with respect to the complex effective field z{£). Without disorder (A = 0), the solution 
is trivial: z' — \ and z" — 0. For finite disorder, the solution can be found numerically in the general case and 
analytically in some limiting cases. 



A. General case 

The numerical solution of Eq. (|l^) for the real and imaginary parts of the effective field for different values of 
disorder in the whole energy range is presented in Fig. |3[ As seen from Fig. ||(a), the real part of the effective spring 
constant varies around its crystalline value z' = 1, being less than unity in the lower part of the energy band and 
greater than unity in the upper part of the band. This is the expected behaviour, because the value of z' de scribe s 
the level-repelling effect for the bare crystalline states when disorder is introduced in the systemn3 (see Sec. [V Ij| ). 
The deviation of z' from unity increases with increasing disorder, which reflects the more pronounced degree of level 
repelling in more highly disordered systems. 

The imaginary part of the effective field is non-zero and negative only in the band region. It basically reproduces 
the shape of the VDOS (Eq. (||)) and describes the widths of the spectral- densi ty peak (Eq. (^) which is proportional 



to the inverse life-time of plane waves in disordered structures (see Sec. IVC). The magnitude of z" increases with 
increasing disorder, thus indicating the shortening of the plane-wave life-time in strongly disordered lattices. The 
imaginary part of the effective field approaches zero in the low-energy limit, z" {e — > 0) — > 0, only if the disorder is 
low enough, for A < A.^, when the structure in equilibrium is mechanically stable within the CPA (i.e. there are no 
negative eigenvalues). 



B. Low-energy limit 

Using the energy dependence of the effective field in the low-energy regime (see Appendix B), the expression (||) 
for the VDOS can be rewritten as 

1 cryst ( ^ \ ^ 1 r,^ryst/^N ^ XDeb „(D/2)-l 



^ - z'{e)' \z'{e)) - i^'iO)f/'' - (z'(0))^/^ ' ^ ^ 



As follows from Eq. (|1^), the low-energy VDOS in systems with lattice disorder functionally behaves the same as in 
the reference crystal, i.e. g"^^^{£) oc e(^/2)-i (see Fig. 0), and differs from the crystalline VDOS just by the factor 



6 



(z'(0)) > 1. In other words, extra states-.iu additiaa to the crystaUine ones, appear in the low-energy regime 
due to disorder-induced level-repeUing efFectsE3L3E3Ejl£3'Lll. The relative density of these extra states is: 

^^^)-^""'(^)^[z'(0)]-^/^-l. (13) 



With increasing energy, the disordered VDOS increases in energy faster than e(^/2)-i and this results in the boson 
peak, the origin of which has been discussed in Rcf.E3. 

In the low-disorder limit (see Appendix B), A ^ 0, the extra density of states (see Eq. (p^) in the low-energy 
range scales quadratically with disorder. 



g"-y«t(e) z 

which has been confirmed numerically (see Fig. H). 



, for A < 1 , (14) 



C. Critical disorder 



All the derivations presented above are valid only for sufficiently small disorder, A < A*. This restriction comes 
from Eq. ( |4ll ) (see Appendix B), which has solutions only for A < A* (see Fig. ^). The critical value of disorder can 
be found numerically, giving A* ~ 1.296. The corresponding value of the real part of the effective field at criticality, 
z'^ — z'(A*), can be found analytically using the fact that {df /dz')A=At = ( see Eq. (pi])), so that 



Eq. (|l|) gives ~ 0.43 for the f.c.c. lattice (with 7 = {Z/2D) - 1 = 1 for Z = 12 and £> = 3). 

When the disorder approaches the critical value where the VDOS behaves as a renormalized Debye function (see 
Eq. (|lj)) and the imaginary part of the effective force constant behaves as z" oc g-^/^ (see Eq. (^)), the low-energy 
range shrinks to zero (see Fig. Above the critical disorder, both the imaginary part of the effective field and the 
effective VDOS are finite at zero energy, z"(0) 7^ and g{e = 0) 7^ 0, and negative eigenvalues appear within the CPA. 
This means that the system becomes vibrationally unstable in equilibrium, even within the mean-field description. 
Such a situation corresponds to the spectrum of the dynamical matrix for instantaneous metastable configurations in 
liquid and glassy states. It turns out that the energy spectrum of the instantaneous dynamical matrix in the unstable 
regime (A > A*) exhibits a peculiar (singular) universal behaviour around zero energy, which we have also found in 
numerical experiments on liquid and vitreous Si02 and the IC-liquid and IC-glas^. 



IV. SPECTRAL DENSITY 



A. Shape of the spectral density 



The spectral density in the plane- wave representation, Akp^s), provides information about the distribution of 
contributions of the different disordered eigenstates in a particular plane wave from the branch /? and with wavevector 
k. In a crystal, the spectral density is obviously a (5-function, A'^p^^is) = S{e — Skp). Disorder introduced in the 
system transforms A'^p'^^i^) into the disordered spectral density which has a shape given by Eq. (^). The shape of 
^k/3(£) depends on the energy behaviour of the effective field z{e). If we assume that the effective field is energy 
independent, then the spectral density has the shape of a Lorentzian peak, the width of which is proportional to 
the imaginary part, z" , of the effective field. Of course, the effective field does depend on energy (see Fig. |3|), but 
if the imaginary part of the effective field is small enough, which is true in the whole energy range for sufficiently 
small disorder, and at least in the low-frequency regime (sk^a ^ see below) at any disorder below the critical one 
(A < A*), then the disordered spectral density still has the shape of a narrow peak (see Fig. |^). The position, e^^^, 
and full width at half-maximum. Ft/?, of such a spectral-density peak can be estimated from the following equations: 

=^ ^'(^kD^k/s and r^p^2\z"{e^^-)\eu0 . (16) 
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For a weak energy dependence of the real part of the effective field (see Fig. ||(a)), we can approximate Eq. ( [l^ ) 
giving: 

e^^- ^ z'(ek^)ek/3 • (17) 

In the low-energy regime (ek/3 0), the imaginary part of the effective field is z" oc e^/^ according to Eq. (j^), so 
that 

-p {D/2) + l 

Tk/s OC . (18) 

This equation corresponds to the well-known Rayleigh law for the peak width, Ti^{uj]ip), in the w-representationBEllEil: 

r^(^k;3) = ^rk/3(e^^) (X u;^;' , (19) 

or r^(a;k,3) c>c to^p in the three-dimensional case (see Fig. ||). 

Eq. (|l^) defines the dispersion law in disordered lattices, lo^^^ = ^/s^p^ — c'^^k (for 0), which is basically 

the same as in the reference crystal except for the renormalized sound velocity, ~ z' (0)^^^^ , with c^'^'*"'* being 
the sound velocity for branch /3 in the reference crystal. Bearing in mind that z'{0) < 1 (Fig. ||(a)), it is evident that 
< c^p^'^^ (see the inset in Fig. |). We have found a similar relation for the model of the topologically disordered 
IC-glass, where c'^^^/c"^^^ ~ 0.72E3. This value corresponds to A ~ 1.28 (see the inset in Fig. ^), for which the 
distribution of the diagonal elements in the diagonal blocks for the disordered-lattice model fits very well the similar 
distribution for the IC-glass (see Fig. 0(a)). 



B. Level-repelling effects 

As seen from Eq. (0), disorder results in a shift of the crystalline spectral density, characterized by the real part 
of the effective field, z'. The bare level moves downwards (upwards) if z' < 1 (z' > 1), i.e. when it is located in 
the lower (upper) part rOfr-t|lie|rbAnd (see Fig. ||(a)). The shift of the bare level, ek/3, by disorder is a consequence 
of level-repelling effectala'Ea'HE^. Indeed, disorder introduced in the system results in the appearance of non-zero 
intera .cJ . iQ n matrix elements between the bare states. These interactions lead to standard quantum repelling between 
levelsEjCa. If the interaction matrix elements are approximately the same for all levels, then the bare crystalline level, 
£k/3, from the lower part of the energy band is shifted downwards just because the number of the states above the 
bare level (and repelling it downwards) is larger than the number of states below it (and repelling it upwards). This 
explains qualitatively why the value of the real part of the effective field is z' < 1 in the lower part of the band and 
vice versa, z' > 1 in the upper half of the energy band (see Fig. |^(a)). 

We have proved quantitatively (analyticallv-and numerically) such a picture of level-repellmg effects in the VDOS 
for the force-constant-disordered fee. latticeEj and iii-^he topologically disordered IC-glassEiJ. A similar effect has 
also been found numerically in models of disordered SiE3. 

An example of the spectral densities for different energies of external plane waves at relatively high disorder (A = 1) 
is shown in Fig. |^. Their widths at sufficiently, large ek/3 become comparable to their positions and plane waves at 
such energies are very short-lived quasiparticlesoE3. The level-repelling effects are clearly seen in pronounced shifts 
(to lower energies) of the peak positions with respect to ek/3. 

We have also presented in Fig. ^ the-cesults of precise numerical solutions for the spectral densities obtained by the 
the kernel polynomial method (KPM)lII. The very good agreement between the CPA and KPM results supports the 
reliability of the mean-field treatment in obtaining the spectral densities. 



C. loffe-Regel crossover and localization 

Disorder also broadens a bare level, so that the 5-functions making up the crystalline VDOS are broadened into 
peaks of finite width, related to the imaginary part of the effective field (see Eq. (|l6|)). This occurs because a plane 
wave |k, /3) is not an eigenstate in the disordered lattice. So it decays with time and can be treated as a quasiparticle 
having a finite life-time. This life-time is invepaely proportional to the width of the peak of the spectral density and 
this determines the physical meaning of z"{e)c3. 

The scattering properties of external plane waves in disordered lattices strongly depend on their bare energy. If the 
energy is low enough, the spectral-density peak is narrow compared to its position, and the quasiparticles live for a 
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long time, thus indicating a weak-scattering regime. For higher energies, the width of the spectral density becomes 
comparable to the peak position (see Fig. |^). This means that the quasiparticles there decay quickly and the bare 
plane waves are in the region pf-slxong scattering. The crossover between weak- and strong-scattering regimes is 



called the lofFe-Regel crossovertBEJS. The energy of the bare plane wave, e]3^, at which such a crossover occurs is 



the lofFc-Regel crossover energy. _ 
The quantitative condition for the lofFe-Regel crossover iso, 



-^r^P - 1 , (20) 



with Lu^^^ = y^£^3 standing for the frequency and Tk/3 being the life-time of the quasiparticle. The quasipar- 
ticle life-time can be easily evaluated in the weak-scattering regime from the time dependence of the probabil- 
ity, I (k, /3; t|k, /3; 0) p, to find the system at time t in the same plane-wave state |k, /?) as it was at t — CH, i.e. 
(|(k,/3;i|k,/3; 0)|2) ~ cos^KTi) exp(-VTk/3), where 



Tk,3 - 777777— Y7^/ ——— oc e^^ . (21) 



The inverse life-time, r^^p, of course, coincides with the full width, Ti^{ekf3), of the spectral density in the uj- 
representationEl, t^^ — T^lskp) = rkp/{2y/e^p^), and the loffe-Regel criterion ( po|) can therefore be rewritten 



[jreseuLHLioiJi_i, r ^ 

as 



2^ r„(Ek;3) . (22) 

This equation should be solved with respect to e^p (or uj\^p = y^Ek/s) and an estimate, (or w^^), for the loffe-Regel 
crossover energy (or frequency) can be found. Such a solution can exist at sufficiently large Wk/3- Indeed, in the low- 
frequency limit, the width of the spectral density, F^^ oc w^^^, is much smaller than the peak position, uj^p^ oc tok/i, 
but F^ grows quickly with increasing Wk/S and can easily reach the peak-position value, if the disorder is large enough 
(see Fig. I). 

We have solved Eq. ( P2|) numerically using the exact energy dependence for the effective field for different values 
of disorder, A (see FigT^). The results for the loffe-Regel crossover energy are presented in the phase diagram 
for vector vibrations in the f.c.c. lattice (see Fig. ||). There is a minimum value of disorder (Amin — 0.7 for the 
disordered f.c.c. lattice using Eq. (p2[)), starting from which the region of the strong-scattering regime appears in the 
system. The strong-scattering regime becomes broader with increasing disorder, while the low-energy weak-scattering 
regime shrinks basically to zero around the critical disorder. A*. We should stress that the boundary between weak- 
and strong-scattering regimes should really be thought of as a relatively wide crossover strip, and Eq. (|2^) serves 
for an order-of-magnitude estimate of the location of this strip. Bearing this in mind, the high-frequency region at 
sufficiently large disorder should not properly be regarded as a region of weak-scattering because the width of the 
spectral-density peak there is of the same order of magnitude as the peak position. The true weak-scattering regime 
occurs only either in the low-energy regime at any disorder less than the critical value, A < A.^, or for all energies 
at small disorder (A < Amin). The existence in these regions of small values of the energy and disorder parameters 
forces the spectral-density width to be much smaller than the peak position. In order to illustrate how sensitive 



the loffe-Regel crossover energy is to the numerical coefficients in Eq. (22), we have calculated from Eq. (|2^) in 
which the spectral-density peak width, F^^, has been replaced, say, by 2Fj^, thus enhancing the strong scattering. The 
results for this loffe-Regel crossover boundary are shown by the dot-dashed curve in Fig. ^. It is clearly seen how the 
region of the strong-scattering regime increases, with just the low-energy and low-disorder regions surviving for the 
weak-scattering regime. i— . 

We have also shown in Fig. ^ the trajectory of the boson peak (circled solid line) calculated as in Ref.B Remarkably, 
it almost coincides with the trajectory of the loffe-Regel crossover for large enough disorder pKhen the strong-scattering 
regime appears in the system. This feature has also been found in real glasses, e.g. v-Si02E3 and in the topologically 
disordered IC-glass model (to be published elsewhere). 

The solid line in the phase diagram shows the upper band edge of the VDOS calculated within the CPA. This appears 
to be a very good estimate for the trajectory of the threshold for the localizaiion-delocalization transition, which, in 
fact, we have calculatecS more exactly by means of a multi-fractal analysia£3 with a relatively high precision, < 5% 
to be situated at slightly lower energies than the CPA band edge. Tie CPA fails to reproduce the high-energy tail in 
the VDOS (see Fig. ||) and more accurate methods should be usedo'EI. It should also be mentioned that, for 1, a 
low-energy tail of localized states extending into the negative eigenvalue range appears in the true spectrumEj. The 
CPA is not able to reproduce this effect for 1 < A < A* but still, for the majority of states with energies far enough 
away from both localization thresholds, the CPA results are very reliable. 
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D. Spectral density for large momentum transfer 



In the previous subsections, we have analysed the situation at small and intermediate energies of external jplaae 
waves. These energies correspond to a relatively small momentum transfer in scattepiig,|experiments (see e.g.oEj). 
How^j,:acy,l4ujgef-|i£a of k are also of interest both from experimental (see e.g.EBCllH'E^) and theoretical (see 
e.g.ll3c3l^-E2ll™lL30) points of view. Attention has been mainly focused on the behaviour of the dynamic structure 
factor, Sjs^ftho), which can be roughly approximated in the high-temperature regime {fuo^/T ^ 1) by the following 
expressionE^: 

S^,{,^)^—-Y,\{dm?5{u^u:^) , (23) 

d 

where the Dcbye- Waller factor has been ignored (cfEi). In Eq. (Eq), we have also neglected the mixture of polarizations 
for phonons |k/3) and assume that /3 is related to the longitudinal branch. In this case, the dynamical structure factor 
is obviously (see Eq. (Q)) proportional to the spectral density, 

where is the spectral density in the a;-representation. 

A peculiar feature has been found in the frequency dependence of the dynamical structure factor in topologically 
disordered glasses at large values of the momentum transfer. Namely., tho .. ap pearance of a secondary peak below 
the Brillouin peak in the w-deoeiidence of S'k/3(w) has been observedtaSBO with increasing value of k. This has 
been related to the boson pe|akO, or interpreted as a manifistation of a microscopic relaxation process typical of a 
topologically disordered glassed. 

A natural question is whether or not the secondary peak exists for the disordered-lattice models and, if so, then 
what is the origin of this peak. In order to answer these questions, first, we plot in Fig. ^(a) the evolution of the 
spectral density Ak;3(e) for the longitudinal branch with increasing value of the wavevector, taken for definitness to be 
along the [100] direction. It is clearly seen from this figure that the spectral density has a pronounced peak (Brillouin 
peak) for all values of k taken from the first Brillouin zone. However this peak is no longer Lorentzian for large 
k. It is asymmetric and has a pronounced shoulder in the low-energy range. A similar picture has been observed in 
simulations for topologically disordered systems. The Brillouin peak shifts to higher energies, becomes vepz, brQad_and 
a smooth shoulder-like feature appears below this peak and may also contain broad peak-shaped featureg£j'Ej'00, so 
that the spectral density starts to resemble the VDOS. In the low-energy regime, the spectral density goes to zero. 

The dynamical structure factor differs from the spectral density and contains an additional decaying energy 
(frequency) function (see Eq. (p^. Therefore, if we plot the effective dynamic structure factor, ^k/3(w) = 

{2mh/T)S)il3{bj) — {k"^ /uj'^)A'^p — [k^ /2^/e)A]^p{e), versus frequency (see Fig. [l^(b)), the low-frequency features 

are enhanced due to the e~^/^ contribution. The low-energy shoulder in the spectral density of the disordered-lattice 
models transforms, for large k values close to the Brillouin-zone boundary, to the secondary peak in the dynamic 
structure factor (see the solid line in Fig. |l^(b)). The intensity of the secondary peak depends on the degree of 
disorder, so that it becomes very pronounced around the critical disorder (see Fig. |ll|(b)), although it does not appear 
for the spectral density in the energy representation (cf Fig. pT|(a)). Therefore, the presence of the secondary peak in 
<S'k/3(w) of disordered-lattice models for large values of wavevector and large degrees of disorder is a distinctive feature 
of such models. Consequently, the occurrence of-a similar peak in topologically disordered structures can hardly be 
attributed solely to topological disorder (cf. RefJl3). 

The origin of the secondary peak in the dynamic structure factor of disordered-lattice models can be understood 
from the qualitative analysis of Eqs. (^) and (]2^). In the low- frequency limit, w — > (or e ^ 0), the energy dependence 
of the spectral density in Eq. (^) is dictated by the imaginary part of the effective field, ^k/3 fx |2"(e)| oc e^/'^ — lo'^ , 
so that S'k/? oc Ak^/cj oc w^. The low-energy limiting behaviour of z"[e) deviates approximately around the energy 
corresponding to the shifted (due to the level-repelling effect) lowest van Hove singularity, i.e. the first minimum in 
the energy dependence of ^'(e) (cf. Figs. |^(a) and (b)). Above this energy, the energy dependence of the imaginary 
part of the effective field is appreciably weaker and basically reproduces the shape of the VDOS (see Fig. H(b)). This 
results in a weaker dependence of the dynamic structure factor, which shows, in this intermediate frequency range, 
either the shoulder-like behaviour, if the wavevector magnitude is not large enough and the low-frequency tail of the 
Brillouin peak plays a significant role (see the long and short dashed lines in Fig. |l^), or a secondary peak for large 
values of k (see the solid line in Fig. |l^), when the influence of the Brillouin peak is suppressed by the function 

~ ' " ' ~ " ' - - ~ - - approximately 

1 However, we 



in Eq. (24). Therefore, the position of the secondary peak in S'k/3(w) for the disordered lattice mode 



coincides with the shifted lowest van Hove singularity and thus with the position of the boson peal 
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should stress that the presence of the pronounced secondary peak in the dynamic structure factor for large values of 
k does not mean a major contribution of plane waves with such k to the boson peak, which is proportional to the 
sum of all plane waves (see Ref.E2l for more detail). 

The cj^-dependence of the dynamic structure factor in the low-frequency regime found for the disordered-lattice 
models is consistent with a /c^-dependence of the spectral, liensit.y width and is in contrast to the constant (non-zero) , 
limit of S'k^(t^) found for topologically disordered glasseal3c3i£5l, and thus can be considered as a drawback of lattice 
models. 



E. Sum-rule correlations and plane-wave scattering 



We have seen in the previous subsection that, independent of the degree of disorder, low-energy external plane 
waves are weakly scattered in a force-constant disordered lattice. This is a consequence of the essentially different 
dependences on the plane- wave energy of the peak width (F^^ oc e^^^ ) and peak position {e^l^^ oc ek/3) of the spectral 
density (see Eqs. ( p7| , p^ ), which gives rise to the occurrence of weak scattering (Ft/? <C ej^f*^) for ek/3 0. The 
question is how generic is this feature. 

In order to answer this question, let us consider a more general formulation of the problem using the following 
Hamiltonian in the site representation: 

H = ^(£. + 7 5]«.,)N)(*|- ^ . (25) 

The diagonal elements of the Hamiltonian have two constituents: random on-site "energies", Ei, and the correlation 
contribution proportional to the sum of random off-diagonal elements, k^ , in the same row. The parameter 7 controls 
the sum-rule correlations between diagonal and off-diagonal elements. If 7 = 1 and Ei — 0, the Hamiltonian ( [25| ) 
describes the scalar vibrational problem, while tke case 7 = corresponds to the well-known electron Anderson 
Hamiltonian with on-site and off-diagonal disordeiH. The spectral properties of a similar Hamiltopian, but defined on 
a positionally random set of sites with deterministic transfer integrals, have been studied in Ref.c3. 

The scattering properties of external plane waves in the disordered structure described by the Hamiltonian ( [25| ) 
can be studied using spectral densities for which the mean-field approach works very well (as shown above for the 
vibrational problem). It is easy to show that only in the case of 7 = 1 and Ei = can the Hamiltonian (^5|) be 
reduced to the form of Eq. (||) which has a multiplicative effective mean field (the energy of the quasiparticles is the 
product of the effective field and the bare energy, ek/3 = 2^(£)ek/3 - see Eq. (p2)). For 7 7^ 1, the situation is quite 
different, because now the effective field has an additive constituent as well, and this drastically changes the scattering 
properties of plane waves. Below, we demonstrate this, considering the simplest particular case for the Hamiltonian 
( p5| ) characterized by 7 = 0, Kij = kq and random Ei taken from the box distribution of half- width, A°', centred 
around Eq = 0, i.e. the standard Anderson Hamiltonian with diagonal disorderEj, for the simple cubic lattice. 

The CPA spectral density, A^{e) = — (7r)~-'^Im(£ — ek)~^, for the Anderson Hamiltonian can be expressed via 
the quasiparticle energy, ek = £k + which is now the sum of the crystalline energy, ek, and the effective field, 
£{e) = e'(£)-fie"(e). The spectral densities A'^{e) are peak shaped and in the low-energy hmit, Se^ = (sk — emfn*) ~^ 
(the energy fck is referred to the lower crystalline band edge, e^f^' = —kqZ, where Z is the coordination number), 
the position and width of the peak can be estimated from the following relations: 



~ fek , (26) 



and 



F^H^k) ^ i7rx'='(A'=')2y&^(x (A<=i)2y&: , (27) 
where the peak position, (5e°',ax = (Smax ~ ^min"^)' measured from the CPA lower band boundary, ej^fj"'^ (see Fig. |2]). 



The last expression is a consequence of the fact that Im [G"y^*{E — e)] ~ —T^X^^y ^ ~ ^mfn*' ~ ^mfif) ~^ ^ with 
being the constant coefficient similar to the coefficient XDcb in Eq. (^7|). 
In full analogy to vibrations (see Eq. (^2|)), the loffe-Regel crossover between weak- and strong-scattering regimes 
for electrons can be defined in terms of the following condition: 

femax^(gk) _ Y''\Ek) . (28) 
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In contrast to vibrations, the spectral-density peak width for electrons scales more slowly than the peak position as 
a function of the plane- wave energy (see Eqs. (2^)-(|2^) and Fig. 12). This means that, in the long- wavelength limit, 
Ssk ^ 0, the width of the spectral density is always larger than the peak position, and the strong-scattering regime 
occurs for plane waves having energies around the lower crystalline band edge (see Fig. |l3|) . The situation is opposite 

for vibrations, for which Fk/j oc e^^^ (cf. Eqs. (p^ ) and (p7|)). As a consequence, the shapes of the regions associated 
with the weak- and strong^cattering regimes for plane waves in the phase diagram are different for electrons and 
vibrations (cf . Figs. ^ and [l3| ) . Now, for electron plane waves with energies around the lower crystalline band edge, the 
scattering is always strong and the weak-scattering regime is realized only in the midband range for sufficiently small 
disorder. In contrast to vibrations, only a single parameter, the strength of disorder, A°', rather than two parameters, 
the degree of disorder A and the energy, for vibrations, controls the appearance of the weak-scattering regime. 
The boundary between the weak- and strong-scattering regimes (see the dashed line in Fig. |l^) is again not precise 
and should be thought of as an estimate for the crossover energy. 

It should be mentioned that the electron phase diagram presented in Fig. O, for the case of very small disor- 



der (A''' 0), supports a standard concept for electronsO, according to which the loffe-Regel crossover from the 
weak- to the strong-scattering regime corresponds to the transition from extended states to localized ones. Indeed, 
the loffe-Regel crossover boundaey (the dashed line in Fig. ^3|) almost coincides with the localization-delocalization 
threshold (the circles in Fig. ^3|)l3 and the CPA band edge (solid line) as A*'' — > 0. At high disorder, the localization- 
delocalization transition, however, is not associated with the loffe-Regel crossover, because all plane waves for such 
disorder are strongly scattered. 

In the more general case of diagonal and off-diagonal disorder coexisting in the Anderson Hamiltonian (but with 
7 = 0), the most successful homomorphic cluster CPAEZflla results in an effective field which necessarily has an additive 
constituent and the above conclusions still hold. In the case of < 7 < 1, the mean field can be constructed as well 
but it also has an additive constituent which disappears only for 7 = 1 (to be published elsewhere). Thus, we can 
conclude that the existence of the low-energy weak-scattering regime is characteristic only of the vibrational problem 
and it results from the exact sum-rule correlations in the corresponding Hamiltonian. 



V. CONCLUSIONS 



We have analytically investigated classical harmonic atomic vibrations in disordered lattices. The f.c.c. atomic 
lattice with force-constant disorder described by a box probability distribution has been chosen for analysis. The 
distributions of matrix elements in diagonal and off-diagonal blocks in the dynamical matrix are rather similar to 
those of topologically disordered models. The atomic vibrations were treated as vectors and analysed within the 
framework of the single-bond coherent potential approximation (CPA); the results are essentially the same as those 
found from a precise numerical study. The CPA self-consistency equation has been derived for vector models and 
solved analytically in the limiting cases of small disorder and low energies. 

The CPA solution of the problem has allowed us to investigate the spectral properties of the model in terms of 
spectral densities. We have shown that the Rayleigh law for the spectral-density width is intrinsic for vector vibrations 
in disordered lattices. 

We have also investigated the loffe-Regel crossover between weak and strong scattering for vector vibrations. The 
regime of weak scattering takes place at all energies for small degrees of disorder, but in the low-energy regime occurs 
only for large disorder. The existence of a disorder-independent weak-scattering regime for low-energy external plane 
waves is a consequence of the sum-rule correlations in the dynamical matrix, typical for the vibrational problem 
only. We have demonstrated that the loffe-Regel crossover occurs in the boson-peak region and is not related to the 
localization-delocalization transition for vector vibrations in the disordered lattice studied. In contrast, the loffe-Regel 
crossover for electrons (described by the Anderson model with on-site energy disorder in the simple cubic lattice) exists 
only at low values of disorder, where it is associated with the localization-delocalization transition. At higher values 
of disorder for the Anderson model, the transition is between strongly-scattered extended states and localized states. 

We have found that a secondary peak appears in the dynamic structure factor below the Brillouin peak for sufficiently 
large values of k for the lattice models as previously found in simulations of topologically disordered glasses. Hence 
this feature cannot be due solely to the presence of topological disorder, nor can it be due to relaxations. 



Appendix A 

Here we give the derivation of the self-consistency CPA equation for vector vibrations. It is knowiSS that the 
self-consistency equation for the effective force constant, K(e), in the single-bond approximation can be written as 
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follows: 

/ 5k 



- lu)(ul=0, (29) 
I - 5K{tj\G\ij) / 



where Sk = Kij — k and the effective Green's function, G = 



e-H 



-1 



, is defined for the effective Hamiltonian, 



H = (k/2)^ Our aim is to find the matrix element, {ij\G\ij), for vector vibrations. This can be rewritten 



as 



= 2 ^(f|f )„(fg))„, (g,„,,„, - G,„,,„,) , (30) 

a.a' 

and this can then be expressed in terms of diagonal elements of the Green's function in the site basis for the reference 
crystal. In order to show this, let us first rewrite the matrix elements, Gia,ja' in the plane-wave basis, {|k, /?)}, 



{zk(r|")-rf)} 



1 <:;aPW1VV1, -Ij- )^ 

G,a,ja' = -TT > (n0)a(n/3)a' r , (31) 

where hp is a _D-dimensional unit polarization vector for the branch [3 and the effective energy, ek/j, obeys the relation: 

ek/3 = 2(e)ek/3 = (r|°^n0) (l-exp{ikr^^°^}) , (32) 

m 

with z[e) = k{e)/KQ being the dimensionless effective force constant and the index running over all nearest 
neighbours to the site i, i.e. j — 1, . . . Z, with Z standing for the number of nearest neighbours (all the sites in a 
crystal are equivalent and Eka, in fact, does not depend on i). 
Substituting expression (P) back into Eq. (pOl) we obtain: 



~ ID 

mG\ij) = — {eG{e)-l] , (33) 



Zk 
where 

G(£)^G..(e)^l^G,„,,„(£) = -i-^^^ (34) 

is the average on-diagonal (independent of i) element of the effective Green's function in the site basis. The effective 
Green's function, G{e), can be expressed via the crystalline one, G"^^^{£), as 

G(e) = iG"y"*(e/z) . (35) 
z 

Finally, the expression for (ij|G|ij) reads: 

{ij\h\i2) = [-G""^' (-) - ll , (36) 

ZzKy^l IZ \Z/ J 

and this solves the problem when substituted back into Eq. (|2^). 

The dimensionless effective spring constant z{e) = H/kq fully determines the effective mean field. This is a complex 
quantity and should be found from the solution of the self-consistuaUpauation (p9|). The self-consistent equation ( |29| ) 
for vector vibrations is very similar to that for the scalar modelsE2rLifE3E3. The only difference is the factor D in the 

denominator of Eq. (p9|), with {ij\G\ij) given by Eq. (^6|). For scalar models, this is absent, even for D-dimensional 
lattices with D 1. 

The crystalline Green's function, G"'^^'^{x), with complex [argument x = e/z, can be easily found as an analytical 
continuation from the real axis of its imaginary part (VD0S)E3. After that, the self-consistent equation can be solved 
numerically in the general case, and analytically in some limiting cases (e.g. e ^ and/or A — > 0), and the effective 
field can be found. 

If z{e) is known, the task of evaluating the spectral density becomes straightforward. Indeed, bearing in mind that 
the energy, iufS, of the eigenstates in the effective crystal is connected with the energy, £k/3, of the similar eigenstates 
in the reference crystal according to relation (|3^), it is easy to see that the VDOS and diagonal elements of the 
spectral-density operator in the crystalline plane-wave basis satisfy Eqs. (^)-(^). 
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Appendix B 

Here we show how the self-consistent Eq. (|l^) can be solved analytically in the low-energy limit, e 0, and how 
the effective field, z{e), can be found in this case. In order to do this, first we find the expression for the function 
a(e) = a'{e) + ia"{£) given by Eq. ( pi] ) via z{e) and e. This can be done, bearing in mind the following relations for 
the crystalline Green's function, C^^'^^e/z), in the low-energy limit, when the Debye approximation for the crystalline 
VDOS, 



cryst 



(D-2)/2 



(37) 



is valid: 



Substituting Eq. 



and 



Re 



^cryst 



(J) 



const and Im 



^cryst ^ — ^ 



into Eqs. (|ll|) and (p6|), we obtain the leading terms in e for a: 

'I \ 

a(e)--^, 



a"{e) - 



2D 



3cb / ^ y 



(38) 



(39) 



(40) 



These expressions can then be used for an expansion of the logarithm in Eq. (|To|) in terms of e, which, after some 
algebra, results in the following final relations for z'{£) and z"{e): 



2A 



z'(7 + l) 



-In 



1 + A + 7z'' 
1 - A-}-7z' 



= , 



(41) 



and 



z"{e) 



2t:Dx 



Deb 



Zz' 



A' + ((7 + 1)/)' - (1 + - [(7 + 1)(^') 



(42) 



where 7 = {Z/2D) - 1. 

In fact, relation ( ^ ) is the equation to find the real part of the effective force constant, z'. The important point is 
that all the parameters entering this equation are energy independent, so that z'{e — 0) ~ ^^'(0) — const. Eq. ( |4l| ) 
can be solved numerically in the case of arbitrary disorder and analytically in the small-disorder limit. 

The imaginary part of the effective force constant, as follows from Eq. (|4^), scales with energy as z"{e) oc e^/^. 
Such a behaviour is basically a consequence of the Debye law for the VDOS in the low-energy regime (see Eq. (|3^)) 
and, in this sense, is general for disordered lattices. This dependence also determines the scaling behaviour of the 
spectral-density peak width with energy, the well-known Rayleigh law in this case. 

Eqs. (^)-(^) have important consequences for the energy dependence of the VDOS in the low-energy limit. Indeed, 
in the limit e ^ 0, the expression (||) for the VDOS can be rewritten in the form of Eq. (| l 2 [) . 

In the low-disorder limit, A ^ 0, Eq. ( ^l|) can be solved with respect to z' while Eq. ([42|) can be simplified, so that 



z'{e) ^ 1 



^A^(l + 0(£)), 



(43) 



'(e) _^TO=bD^z,/2^2 (44) 
3Z 



These relations are used in deriving Eq. (p^ 
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FIG. 1: The distribution of the matrix elements (scaled by the spectral band widths) of the vibrational Hamiltonian matrix 
for the IC-glass (solid line) and disordered f.c.c. lattice (dashed (A = 1) and dot-dashed (A = 1.28) lines): (a) p^Kn) for the 
traces of the diagonal blocks; (b) p{Hia,ia') for the off-diagonal elements in the diagonal blocks; (c) p{Kij) for the traces taken 
with opposite sign in the off-diagonal blocks; (d) p{Hia,ja') for the off-diagonal elements in the off-diagonal blocks. The inset 
in (c) shows p{Hf^^^) for nearest-neighbour interactions only. 
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FIG. 2; The vibrational density of states (VDOS) versus energy in the crystalline (dashed line) f.c.c. lattice and in f.c.c. 
lattices with force-constant disorder characterized by different widths (as marked), A, of the box [distribution. The curves for 
the disordered models have been obtained by the CPA and the kernel polynomial method (KPM)lj, as described in the text. 
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FIG. 3: (a) The real, z' , and (b) imaginary, z" , parts of the effective force constant versus energy for different values of 
disorder as marked. 
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FIG. 4: The VDOS, g(e), in the low-energy regime for different values of disorder as marked. The thin solid line (guide for 
the eye) shows an e^l'^ dependence. 
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FIG. 6: The evolution of the function f{z') given by Eq. ([4l|). The larger root of the equation f{z') = existing only at 
A < A* corresponds to the real part of the effective force constant. 
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FIG. 7: The spectral densities for vector vibrations in the f.c.c. lattice with force-constant disorder (the width of the box 
distributions A = 1) calculated by the CPA (solid lines) and KPM (dashed lines). 




FIG. 8: The peak positions, cu'^y^ (thin linos), and peak widths, 2TTr^ (thick lines), versus plane- wave frequency, iOka, for 
different values of disorder as marked. The crossing points of the curves for peak positions and peak widths correspond to the 
loflfe-Regel crossover frequency. The thin solid line (guide for the eye) shows an ui^/^ dependence. The inset shows the ratio of 
the sound velocities in disordered and ordered lattices, cjj'^/cj^''', versus degree of disorder (solid line). The dashed curve in 
the inset corresponds to the low-disorder approximation, cj^"'/cj^''' ~ 1 — DA^ /3Z, valid for A <C 1. The critical behaviour of 
^'74''''* is evident for A ^ A,. 
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FIG. 9: Phase diagram for vector vibrations in the f.c.c. lattice with force-constant disorder characterized by the box 
distribution. The dashed (dot-dashed) line corresponds to the loffe-Regel crossover plane-wave energy obtained from Eq. ( |2^ ) 
(or from Eq. ( p^ ) in which the spectral-density peak width is decreased, say, by a factor 2). The circled solid line shows the 
trajectory of the boson peak (for A = 0, its position coincides with the lowest van Hove singularity in the crystalline lattice). 
The solid line represents the upper band edge obtained by the CPA. The horizontal dotted line separates the systems stable in 
equilibrium from those that are mechanically unstable. 




FIG. 10: The evolution of the spectral density, j4k^(e) (a), and dynamic structure factor Sk0 (b), with wavevector magnitude, 
k, as marked (the corresponding crystalline energies, £k/3 are given in brackets), taken along the [100] symmetry direction in 
the first Brillouin zone (the maximum value of fc = 2.72 corresponds to the zone boundary) for the longitudinal branch in the 
f.c.c. disordered (A = 1) lattice. The thin solid line in (a) for k = 2.72 shows the numerical results obtained by KPM. 
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FIG. 11: The evolution of the spectral density, Ak/3(e) (a), and dynamic structure factor Sul/b (b), with degree of disorder, A, 
as marked for k = 2.72 (longitudinal branch) corresponding to the zone boundary along [100]. 




FIG. 12: The peak positions, (5emax (thin lines), and peak widths, 27rr''' (thick lines) for the electron spectral densities for 
the Anderson Hamiltonian with on-site disorder for the simple cubic lattice, versus plane-wave energy, 5eu, at different values 
of disorder as marked. The crossing points of the curves for peak positions and peak widths correspond to the loflfe-Regel 
crossover energy for electrons. 
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FIG. 13: Phase diagram for electrons described by the Anderson model with on-site energy disorder in the simple cubic lattice. 
The dashed line corresponds to the loffe-Regel crossover plane- wave energy obtained from Eq. (^). The solid line represents 
the lower band edge obtained by the CPA and the circles show the trajectory for the localization-delocalization thresholdEI. 



